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Abstract 

We propose a model for the production and suppression of J/ip mesons in 
high-energy hadronic collisions. We factorized the process into a production 
of the cc pairs of relative momentum k = k c — kc convoluted with a transition 
probability from the produced cc pairs into the observed J/tjj mesons. As the 
produced cc pairs exit the nuclear matter, the multiple scattering between the 
colored c and c and the nuclear medium increases the square of the relative 
momentum between the c and c, q 2 = —k 2 , such that some of the cc pairs 
gain enough invariant mass to transmute into open charm states. With only 
one parameter, the amount of energy gained by the produced cc pair per unit 
length in the nuclear medium, our model can fit all observed J/ip suppression 
data including recent NA50 data from Pb-Pb collisions. 

PACS: 12.38.-t, 12.38.Bx, 13.85.Ni, 14.40.Lb, 25.75.-q 
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The suppression of J /if; production in high energy nucleus- nucleus collisions has been 
suggested as a potential signal for the existence of the quark-gluon plasma [[|]. In recent 
years, strong J /if) suppression has been observed in high energy hadron- nucleus and nucleus- 
nucleus collisions by a number of experiments and various theoretical explanations of 
the observed J/tp> suppression have been proposed 0. Recently, the NA50 Collaboration 
at CERN observed a much stronger J/if) suppression in Pb-Pb collisions at SPS energies 
||. It has been argued that a "conventional" approach cannot explain this new data, and 
debates and controversies on the theorized origin of the suppression have reached a new level 
|J. In this letter, using the same mechanism proposed in Ref. J7| for the J/if> suppression 
in hadron-nucleus collision, we explain the newly observed strong J /if) suppression in high 
energy Pb-Pb collisions. 

The production of J/if> mesons in high energy hadronic collisions is believed to have two 
factorizable stages: the production of cc pairs, and the formation of J /if) mesons from the 
produced cc pairs. Due to the large mass of the charm quark, the cc pair production should 
be a short- distance process in nature, and therefore calculable with perturbative methods. 
On the other hand, the formation of J/if> mesons from the initially compact cc pairs takes 
a relatively long time and is nonperturbative ||. Throughout the years, the main debate 
on the mechanism of J /if) production has focused on the second stage. Three models are 
commonly used in the literature for calculating the cross sections of J /if) production: the 
Color Evaporation Model (CEM) @, the Color-Singlet Model (CSM) @, and the Color- 
Octet Model (COM) |pT| . Since the data on J/ip suppression in nucleus-nucleus collisions are 
presented in terms of total cross sections f4j||, we will focus the remainder of our discussions 
on the total production rate of J /if) mesons in hadronic collisions at fixed target energies. 
For calculating the total J /if) production rate, we propose a simple model which is general 
enough to admit any of the three possible formation mechanisms mentioned above. 

For the collisions between hadrons (or nuclei) A and B, A(pa) + B(p b ) — > J/ip(Pj/^)+X , 
the total J/ip cross section can be factorized as follows, 
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where the four-vector = + k^ is the total momentum of the produced cc pair. In 
Eq. ([I]), the variable q 2 is equal to the square of the relative momentum between the c and 
c in their rest frame, q 2 = (2k c ) 2 . If the c and c can be approximated as on their mass-shell, 
kc = kd = m c> an d we have q 2 = Q 2 — Am 2 . Without trying to separate contributions from 
different color channels, we define the F c5 ^j/^(q 2 ) in Eq. ([[]) to be a transition probability 
for a color-averaged cc pair of the relative momentum square q 2 to evolve into a physical 
J /if) meson. We propose three alternatives for parameterizing the transition probability, 
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where m' is the mass scale for the open charm threshold. In Eq. (0), Nj/^ and are 
to be fixed by fitting the existing total production cross section data from hadron-hadron 
collisions. 

The transition probabilities in Eq. (fj) represent a wide range of J/ip formation mech- 
anisms. The F (c) (g 2 ) implies that all cc pairs with invariant mass below the open charm 
threshold have the same constant (C) probability to become the physical J ftp mesons, which 
is effectively the same as the Color Evaporation Model M. The F^ G \q 2 ) corresponds to fol- 
lowing assumptions: the transition amplitude (cc\J/ip) does not involve any radiation and 
interaction with the medium, and it is then proportional to the J/i/j wave function parame- 
terized as a Gaussian (G). If we neglect the g 2 -dependence in the production of the cc pairs 
m Eq. (H|), and require the cc to be color- singlet, the total cross section with F^ G \q 2 ) is 
effectively the same as that from the Color-Singlet Model [|Kj . 

If the J /ip mesons are formed after a long-time expansion from the small size cc pairs, 
and radiating soft gluons adjusts the color of the pairs, it is then natural to assume that 
the g 2 -dependence of the transition probability is associated with that radiation, and to 
choose a power-law (P) distribution, F^ p \q 2 ) in Eq. (|2c|), for the transition probability. If 
we expand the transition probability at q 2 ~ 0, the normalization of F( p \q 2 ) can be related 
to the combination of the matrix elements in the Color-Octet Model . 

In principle, with a different functional form of F(q 2 ), our factorized formula in Eq. (|]) 
can be generalized to calculate the total cross sections for producing other quarkonium 
states. 

To evaluate the J/ip total cross section in Eq. (|TJ), we need to calculate the production 
rate for the cc pairs at invariant mass Q 2 . As argued in Ref. [Q, the production rate can be 
factorized into a convolution of two parton distributions from the two incoming hadrons and 
a short-distance hard part, da a+ b^ C c+x / dQ 2 , which represents the perturbatively calculable 
hard parts for the parton a and b to produce the cc pairs with mass Q 2 . Similar to the 
total Drell-Yan cross section, the one-scale cross section da/dQ 2 for producing the cc pairs 
at the fixed target energies should be well-represented by the leading order calculations in 
a s , and the high order corrections are given by a smooth K-f actor. At the leading order in 
a s , the partonic contributions come from two subprocesses: qq — > cc and gg — > cc. With the 
K-factor for effective high order contributions, the total J/ip cross section in Eq. (|J) can be 
written as M 



cta+b^j/iP+x(S) = Kj/ip ^2 J dq 
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where Y^ a ,b runs over all parton flavors, and Q 2 = q 2 + Am 2 . Because of the two-parton 
final-state at the leading order, the incoming parton momentum fractions are fixed by the 



kinematics, and given by x a = (yx 2 7 + AQ 2 / S + Xf)/2 and Xb = {Jx% + AQ 2 / S — xp)/2, 
respectively. In Eq. (^[), the short- distance hard parts for producing the cc pairs, a g g^ C c(Q 2 ) 
and cr gg ^ C c(Q 2 ), are given in Refs. ||7|jl3[|. In Eq. (|3]), the integration limits of x F are chosen 
to be consistent with the data, and the limits of q 2 are specified by the functional form of 
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F c5 ^j/^(q 2 ) in Eq. (0). Eq. ([3]) combining the transition probability defined in Eq. (|]) is our 
model for calculating the total J/ip hadronic cross sections. 

In Fig. 0, we plot the total J/i/j cross sections using Eq. (|3|) in comparison with the data 
in hadronic collisions |14[]. Following the same fitting approach used in Ref. f\A \, we fix all 
parameters in Eq. (0) and list them in Table |, in which /j = Kj Nj ^ is defined as an 
overall normalization factor. To obtain the theory curves in Fig. [3], we used CTEQ4L parton 
distributions [1131, and noticed that EMC effect gives a very small modification to the total 
cross sections because of the integration of xp and q 2 j7|. In addition, we set m c = 1.50 GeV 
and m! = 1.869 GeV. Choosing different values for the m c and m' changes the fitting 
parameters in Table | slightly. But, it does not change the quality of the comparison in 
Fig. |l|. All three parameterizations in Eq. @j provide a good fit to the total J/ip cross 
sections from proton-nucleon collisions at fixed target energies. 

In nucleus-nucleus collisions, the produced cc pairs are likely to interact with the nuclear 
medium before they exit. Observed anomalous nuclear enhancement of the momentum 
imbalance in dijet production tells us that a colored parton (quark or gluon) experiences 
multiple scatterings when it passes through the nuclear medium, and the square of the 
relative transverse momentum between two-jets increases in proportion to the size of the 
nucleus ]TB| , P^7| . If we let the c and c be the parent-quarks of two jets, the q 2 becomes 
the square of the relative momentum between the two jets in their cm. frame. Therefore, 
as the c and c pass through nuclear matter, just like a di-jet system, the square of the 
relative momentum q 2 increases. As a result, some of the cc pairs might gain enough relative 
momentum square q 2 to be pushed over the threshold to become open charm mesons, and 
consequently, the cross sections for J/ip production are reduced in comparison with nucleon- 
nucleon collisions. 

If the formation length for the J /if) meson, which depends on the momenta of the cc 
pairs produced in the hard collision, is longer than the size of the nuclear medium, it is 
reasonable to assume that the transition probability F cE ^j/^,(g 2 ), defined in Eq. (|2|), can be 
factorized from the multiple scattering. Then, as far as the total cross section is concerned, 
the net effect of the multiple scattering of the cc pairs can be represented by a shift of q 2 in 
the transition probability, 



Q 



q 2 + e 2 L(A,B) 



(4) 



In Eq. (IJ), L(A, B) is the effective length of nuclear medium for the cc pair to pass through 
in the collisions of two nuclei of A and B, and it depends on the details of the nuclear 
density distributions [I5|. In Eq. ([|), the e 2 represents the square of the relative momentum 
received by the cc pairs per unit length of the nuclear medium. The value of the e 2 can be 
estimated from the observed nuclear enhancement in the momentum imbalance of two-jets in 

Using the data from pion-nucleus collisions []T3, we estimate 
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hadron-nucleus collisions 
that e 2 ~ 0.2 — 0.5 GeV 2 per unit length of nuclear medium [13,17|. 

In Fig. 0, we plot the predictions of J/ip total cross sections in proton-nucleon, proton- 
nucleus and nucleus-nucleus collisions. The data in Fig. are from Ref. [|19|], in which all 
data were rescaled to Pbeam = 200 GeV. The effective length L(A, B) were chosen to be the 
same as those used in Ref. [19|. Three theory curves correspond to three parameterizations 
defined in Eq. (^). The values of the parameter e 2 for different parameterizations are listed in 
Table |, which are consistent with our earlier estimates from the momentum imbalance of a 
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di-jet system. As in Fig. |], we used the CTEQ4L parton distributions and set m c = 1.50 GeV 
and m' = 1.869 GeV for plotting Fig. |[ Choosing different values for the m c and m' changes 
the values of e 2 slightly, but, it does not change the features of Figs. [I] and 0. 

The three parameterizations of the transition probability in Eq. (FJ) represent the different 
J /if) production mechanisms, and naturally, they predict different behavior of J /if) suppres- 
sion in Fig. Therefore, understanding the suppression can also help us to distinguish the 
production mechanism. 

For the parameterization F^ G \q 2 ) in Eq. (pbj), a shift of q 2 to q 2 in Eq. (|4]) for the J/ip 
suppression in nucleus-nucleus collisions gives following relation 



<?AB-*3/^{S) = exp 



e 2 



2a 2 F 



L(A,B) 



a NN ^j/^(S). (5) 



This relation is effectively the same as that predicted by the Glauber theory, if we let 
the suppression factor in the simple Glauber theory be exp[— cr abs pL(A, B)}, with p being 
the nuclear density. With the parameters in Table |, we have the effective absorption cross 



section cr a b s ~ 5.9 mb, which is the same as that obtained in Ref. [ 19[ . In addition, our model 
interprets the <7 a b s in Glauber theory is proportional to the energy absorbed by the colored 
cc pairs. As expected |J, like the Glauber theory, the parameterization of F^ G \q 2 ) does 
not generate enough suppression for heavy nucleus-nucleus collisions. As discussed earlier, 
the Gaussian parameterization corresponds to assuming that the formation process does not 
have radiation, and is then proportional to the square of the J /if) wave-function. However, 
as we learned from the success of the Color-Octet Model |TTJ, the formation of J /if) should 
involve the expansion of the cc pairs and the radiation of soft gluons. Therefore, we should 
prefer a power-law parameterization than a Gaussian parameterization, and consequently, 
we expect to have more suppression than that expected from the Glauber approach. 

Since the F^ c \q 2 ) is the same as the F^ F \q 2 ) when cuf = 0, we expect a maximum 
suppression from the Color-Evaporation Model, as seen in Fig. [2|. The Color- Evaporation 
Model assumes that all cc pairs with invariant mass less than the open charm threshold 
should have the same transition probability to become the J /if) meson. However, the phase 
space cutoff for the cc pairs to become the open charm in Eq. (^) appears classical. In 
quantum theory, the cc pairs with invariant mass less than the Am' 2 should have a small, 
but non-zero, probability to become open charm systems, and in general, the pairs just below 
the threshold should have a relatively larger probability for a transition to open charm than 
those far below the threshold. With these considerations, we believe that the power-law 
parameterization F^\q 2 ) with a > represents the more accurate physics for the J /if) 
production. 

Our model of the J /if) production relies on the factorization between the production of 
the cc pairs and the formation of the J /if) mesons, and accurate calculation for the production 
of the cc pairs. For the total cross sections at fixed target energies, we believe that such 
factorization is justified, and the perturbative calculations are reliable. However, at collider 
energies, most 3/ip events are measured at large transverse momentum, Qt', and the 
as well as the Qt are much larger than the invariant mass of the cc pairs, Q, in Eq. @. 
Therefore, the perturbative calculations of da/dQ 2 and da/dQ 2 dQ\ for the production of 
the cc pairs become very nontrivial due to the large logarithms log(l/x) ~ log(S/Q 2 ) and 
\og{Q\ / Q 2 ) . A resummation of such large logarithms to all order in a s is necessary in order 
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to have a reliable prediction. We defer our detailed discussions on J /if; production and 
suppression at the collider energies to another publication. 

Our predictions of the J/ip suppression, as shown in Fig. ^], depend on an additional 
assumption: the separation of the multiple scattering of the cc pairs and the formation of 
the J ftp mesons. We believe that this additional assumption can only be justified when the 
J /if) formation length is larger than the effective medium length L(A, B) in our Eq. (f|). Once 
the J /if) meson is formed, the multiple scattering with nuclear medium should be reduced due 
to the color singlet nature of the meson, and then, the Glauber formalism for the suppression 
should be more relevant. Therefore, if there is no QCD phase transition to the quark-gluon 
plasma, we expect the following features for the J/ip suppression in nucleus-nucleus collisions. 
As the size of colliding nuclei increases, the J/ip suppression should follow the dotted curves 
in Fig. |2|; and when the L(A, B) is compatible to the J /if) formation length, the suppression 
will become smaller than what is predicted by the dotted curve. Finally, we conclude that 
our simple model for the total cross sections of J /if) production in nucleus- nucleus collisions, 
as defined in Eq. (||), can explain the existing data in hadron-hadron, hadron- nucleus and 
nucleus-nucleus collisions. Although our model is different in many aspects from what have 
been discussed in the literature, we believe that it has the key physical mechanisms for the 
J/if> production and suppression. Detailed comparison between our work and others in the 



literature are given in Ref. [13 
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FIGURES 




FIG. 1. Total J/tp cross sections in proton- nucleon collisions as a function of the colliding 
energy \/~S. 
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FIG. 2. Total J/ip cross sections with the branching ratio to fi + /j> in proton-nucleon, 
hadron-nucleus and nucleus-nucleus collisions as a function of the effective nuclear medium length 
L(A,B). 
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TABLES 



TABLE I. Values of parameters used to produce the theory curves in Figs. |l] and ^. 





i?(C) 




i?(P) 


fj/lp 


0.248 


0.470 


0.485 







1.2 GeV 


1.0 


e 2 (GeV 2 /fm) 


0.45 


0.29 


0.25 
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